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Abstract. Intrinsic stellar variability can hinder the detection of shallow transits, particularly 
in space-based data. Therefore, this variability has to be filtered out before running the transit 
QQ , search. Unfortunately, filtering out the low frequency signal of the stellar variability also modifies 
• the transit shape. This results in errors in the measured transit depth and duration used to derive 
■ the planet radius, and orbital inclination. We present an evaluation of the magnitude of this 
04 , effect based on 20 simulated hght curves from the CoRoT blind exercise 2 (BT2). We then 
(— ( ■ present an iterative filter which uses the strictly periodic nature of the transits to separate them 
' from other forms of variability, so as to recover the original transit shape before deriving the 
planet parameters. On average with this filter, we improve the estimation of the transit depth 
QQ and duration by 15% and 10% respectively. 
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O ' !• Introduction 

' l.T Planet parameters from observations and associated errors 

Ci , The radius (i?p, eq. II. ip and the mass (Afp, eq. IT2p of an exoplanet can be fully solved 
when measuring both the flux and the radial velocity variations of the parent star due 
^-H ■ to its orbiting planetary companion. 
> • 

'■ [AF 

rn, ^P-^^— (1-1) 
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^ Mv = M* (1.2) 



where, i?^, and are the radius and the mass of the parent star, the flux variation 
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wiiere, n,^, anu iw* are trie lauius aiiu iiic mass ui uiie pareiiu siar, — p- 

• • I due to the planet transiting the disc of its parent star, K the amplitude of the radial 

. ^ velocity variation of the parent star due the gravitational influence of its orbiting planet, 

. i and P the orbital inclination and period of the planet, and G the gravitational constant. 

^ ' P, and i can be measured from the light curve. A common way to measure R^, and 



M^, is by comparing the stellar spectrum to stellar atmosphere models, allowing to derive 
the stellar parameters (Te//, logg, [Fe/H]) used to obtain the stellar mass and radius, i?* 
can also be measured more precisely, and without the use of models, with interferometry, 
or with transit fltting (in the case of high precision light curves). 

1.2. Planet parameters and planet evolution and formation models 

Improving the precision on observational planet masses and radii is important for both 
planet structure and planet formation models. The internal structure of a planet can be 
studied by comparing its mass and radius to model predictions of planets with different 
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composition. Determining planet structure is important to derive observational statistics 
on planet types, which can then be compared to the predictions of planet formation 



models. Seager et al. (2007) show that to determinate the composition of sub-Uranus 
planets, error bars of 2% on the planet parameters are required. The current uncertainties 
on planet masses and radii are of the order of 10%. Improving these measurements is 
thus vital to help confirm the models. 

1.3. Sources of uncertainties on planet parameters 

The uncertainties on the planet mass, radius, and inclination depend on the uncertainties 
on the host star mass and radius, on the uncertainties on the transit parameters (-^-j 
total transit duration), and on the uncertainties on the radial velocity measurements. For 
large planets (> Jupiter), the uncertainties on the planet mass and radius are mainly due 
to the uncertainties on the stellar parameters. For sub-Uranus planets around active stars, 
the uncertainties on the planet mass and radius can be dominated by the uncertainties 
on the transit parameters. 



2. Side effects of pre-detection stellar variability filters 

2.1. Pre-detection stellar variability filter 

Active stars shows intrinsic flux variations due to temporal evolution and rotational mod- 
ulation of structures on their surface (stellar spots, plages, granulation). Intrinsic stellar 
variations can have amplitudes much greater than transits, and thus, can hinder transit 
detection. These variations occur at a lower frequency (longer time length) than tran- 
sit events (chap. 3, fig. 3.2). Pre-detection stellar variability filters use this difference to 
separate variations typical of stellar variability, from those on the time scale of transits 
(minutes to hours). 



transit + instrumental noise 



transit + instrumental noise 





tiltered from stellar variability - Aigrain&lrwin(2004) filter 




Figure 1. Stellar variability hinders transit detection { column cl ): the bottom plot is an example 
of a CoRoT BT2 light curve with Jupiter like transits on an active host star. This light curve 
is composed of a simulated transit signal (top plot), the instrumental noise expected in CoRoT 
data, and simulated stellar variability. Pre-detection filters can deform the shape of transits 
{ column b ): the top plot of column b is the phase fold of the top plot of column a, centered on 
the transit. The bottom plot of column b is the phase fold of the bottom plot of column a, filtered 
from stellar variability using Aigrain & Irwin (2004) filter. 



Stellar variability and planetary transits 
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The work presented in this paper is based on the pre-detection stellar variability filter 



described in Aigrain & Irwin (2004) This filter is a combination of median filtering (to 
reduce the level of random noise), and of boxcar averaging/smoothing (to extract the 
long term variations: the stellar variability). 

2.2. Deformation of the transit signal 



We tested the effect of Aigrain & Irwin (2004) pre-detection filter on the transit sig- 



nal of 20 simulated light-curves, produced for the CoRoT bhnd test exercise 2 (BT2, 



Moutou et al. (2005)). The three main components in these light-curves are the stellar 



variability, the transit signal, and the instrumental noise. Figure[T^ (bottom panel) shows 
a BT2 light curve of an active star with Jupiter-like transits. 

We noticed that 'Aigrain & Irwin (2004) filter deforms the transit shape (Figure [T]d, bot- 
tom panel). This effect is stronger for shallow transits on active stars. On average, for 
filtered light curves, the transit depth is under-estimated by 20% and its duration by 
15% (compared to the value measured from the original transit with instrumental noise 
only) . Miss-estimating the transit depth and duration leads to a miss-estimation of the 
planet radius, and orbital inclination (used to derive the planet mass). 



3. Iterative filtering 

We designed an iterative process to filter stellar variability and extract the original 



transit shape, based on Aigrain & Irwin (2004) filter, and using the additional post 



detection knowledge on the planet orbital period. 
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Figure 2. Top: Chart of the iterative filter process. The light curve (y) is decomposed into three 
components: the transit signal (a), the stellar variability (/), and the residual noise (r). The final 
estimation of the transit signal (oflnai) is obtained after iterating on the pre-detection filter (step 
called "smoothing"). Bottom: Transit (aflnai) resulting from the iterative filter applied to the 
light curve of Figure [1^ bottom plot. 
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3.1. Implementation of the post- detection filter 



Our post-detection filter is based on Kovacs et al. (2005) idea of decomposing the light 



curve into three components (transit signal, correlated noise, and white noise), to filter 
out systematics due to atmospheric fluctuations. In our case, we filter out stellar intrinsic 
variability extracted using an iterative process (top chart of Figure ^ . 



We estimate the stellar variability by applying Aigrain & Irwin (2004) pre-detection fil- 
ter. We subtract the estimation from the original light curve, phase-fold the result, and 
bin it to reduce the random noise. We then remove the resulting transit signal from 
the original light curve, and start a new iteration re-evaluating the stellar variability 
in the light curve freed from the transit signal. We stop iterating when the residuals 
(original light curve minus last estimations of the transits and stellar variability) stop 
evolving. The best estimation of the transit signal is the original light curve minus the 
last estimation of the stellar variability (after 3-4 iterations). 

3.2. Evaluation of performance 
We have applied the iterative stellar variability filter described in section 13.11 on 20 
CoRoT BT2 light-curves with transits. Iterating on the estimation of the stellar vari- 
ability, and using the additional knowledge on the orbital period, appears to successfully 
better recover the original transit shape (compare bottom of Figure [5] with bottom plot 
of Figure [T|d). On average, light curves filtered with the iterative filter gives a better 
estimation of the transit depth by 15% and of the transit duration by 10%, than when 



filtered with Aigrain & Irwin (2004) pre-detection filter. 

The iterative filter has also been found efficient in reconstructing shallow transits that 
had become un-detectable after pre-detection stellar variability filtering. Another use for 
the iterative filter could thus be to confirm borderline transit detections in active stars. 



4. Summary and future work 

Based on 20 simulated CoRoT light curves from the BT2 light curve sample, we have 
shown that Aigrain & Irwin (2004) pre-detection filter, used to remove stellar variability 
prior to transit detection, deforms the transit depth by 20% and the transit duration by 
15%. To circumvent this, we have adapted Kovacs et al. (2005) iterative filtering method 
to the case of filtering stellar variability present in space-based light curves. The resulting 
post-detection iterative filter improves the estimation of the transit depth and duration 
by 15% and 10% respectively. 

The two areas where we plan to focus future efforts are to further automate the filter- 
ing process (some user interaction is currently needed to initialize the filter smoothing 
length), and b) to evaluate the improvement in the planet parameters uncertainties re- 
sulting from the improvement in the transit parameters from this work. As b) depends on 
the particulars of each system, we plan to derive these uncertainties using Monte Carlo 
simulations for known planets. Further tests will include using other pre-detection filters, 
including ones which can be applied to data with significant temporal gaps. 
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